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For hexagonal nets, descriptive of {111} fee surfaces, we derive from combinatoric arguments a 
simple, low-temperature formula for the orientation dependence of the surface step line tension and 
stiffness, as well as the leading correction, based on the Ising model with nearest- neighbor (NN) 
interactions. Our formula agrees well with experimental data for both Ag and Cu{lll} surfaces, 
indicating that NN-interactions alone can account for the data in these cases (in contrast to results 
for Cu{001}). Experimentally significant corollaries of the low-temperature derivation show that 
the step line tension cannot be extracted from the stiffness and that with plausible assumptions 
the low-temperature stiffness should have 6-fold symmetry, in contrast to the 3-fold symmetry of 
the crystal shape. We examine Zia's exact implicit solution in detail, using numerical methods 
for general orientations and deriving many analytic results including explicit solutions in the two 
high-symmetry directions. From these exact results we rederive our simple result and explore subtle 
behavior near close-packed directions. To account for the 3-fold symmetry in a lattice gas model, 
we invoke a novel orientation-dependent trio interaction and examine its consequences. 



PACS numbers: 68.35.Md 05.70.Np 81.10.Aj 65.80.-l-n 
! I. INTRODUCTION 

Most of our current understanding of surface morphol- 
ogy is based on the step-continuum model, ^ which treats 
the step itself as the fundamental unit controlling the 

\ evolution of a surface. In this model, the step stiffness 
/3 serves as one of the three fundamental parameters; it 
gauges the "resistance" of the step to meandering and ul- 
timately accounts for the inertia of the step in the face of 
driving forces. The stiffness can be derived from the step 
line tension (3{9), the excess free energy per length associ- 
ated with a step edge of a specified azimuthal orientation 
0. From (3, the two-dimensional equilibrium crystal shape 
(i.e., the shape of the islands) is determined. 

The goal of this paper is to find low-temperature (T) 
formulas for (3{e) and thence (3{e) = I3{e)+(3"{e) as func- 
tions of the azimuthal misorientation 0, assuming just 
nearest-neighbor interactions in plane and an underly- 

i ing {111} surface. Such surfaces are characterized by a 
six-fold symmetric triangular (hexagonal) lattice, allow- 

\ ing all calculations to be done in the first sextant alone 
(from 0° to 60°). In contrast to I3{9), we shall find that at 

\ low T, j3{0) is insensitive, under plausible assumptions, 
to the symmetry-breaking by the second substrate layer, 
so that plots from 0° to 30° suffice. Although an analytic 
solution exists for the orientation dependence of /3(0) on 
a square latticejSiSi^ only an implicit solution to a 6*'*- 
order equation has been found for a hexagonal latticed 
This makes comparisons to experiment rather arduous, 
particularly when trying to compare data on (3{0), which 
is related to I3{0) through a double derivative with re- 
spect to 9. Fortunately, we will see that a remarkably 
simple formula exists for the orientation dependence of (3 
at temperatures which are low compared to the charac- 



teristic energy of step fluctuations (i.e. the kink energy 
or the energy per length along the step). For noble met- 
als, room temperature lies in this limit, facilitating direct 
comparisons to experiment. 

Motivating this work is a recent finding^ that the 
square-lattice nearest-neighbor (NN) Ising model under- 
estimates /3 by a factor of 4 away from close-packed di- 
rections on Cu{001}. Later work^'®'^ showed that much 
(but not all) of this discrepancy could be understood by 
considering the addition of next-nearest-neighbor (NNN) 
interactions. For the triangular lattice, we will see that 
such a longer-range interaction is not required to describe 
the orientation dependence of (3. 

In the following Section, we characterize steps on a 
hexagonal lattice and perform a low-temperature expan- 
sion of the lattice-gas partition function, assuming only 
NN bonds are relevant, and derive both [3{6) and /3(0). 
We obtain a remarkably simple expression for the latter 
in Eq. H18|l . Since this low-T limit is determined by ge- 
ometric/configurational considerations, it becomes prob- 
lematic near close-packed orientations {0 — Q°), where the 
kinks must be thermally activated. Therefore, we make 
use of exact results to assess in several ways how small 6 
can be before the simple expression becomes unreliable. 
In Section III, we present three general results for island 
stiffness that are valid in the experimentally-relevant low- 
T limit when configurational considerations dominate the 
thermodynamics. We show that the line tension can- 
not be [re] generated from the stiffness and that the stiff- 
ness can have full 6-fold symmetry even though the sub- 
strate and the line tension have just 3-fold symmetry. 
Accounting for such 3-fold symmetry with a lattice-gas 
model on a hexagonal grid requires an extension from the 
conventional parametrization; we posit an orientation- 
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dependent interactions between 3 atoms at the apexes of 
an equilateral triangle with NN legs. In Section IV, we 
compare our results to experiments on Ag and Cujlll} 
surfaces, showing that Eq. (|18|) provides a good approx- 
imation and, thus, that NNN interactions are much less 
important than on Cu{001}. The final section offers a 
concluding discussion. Three appendices give detailed 
calculations of the leading correction of the low-T ex- 
pansion, of explicit analytic and numerical results based 
on Zia's exact implicit solution,^ and of Eq. H18(l as the 
low-T limit of Zia's solution. 



II. ISING EXPANSION ON A TRIANGULAR 
LATTICE 

A. Recap of Results for Square Lattice 

The orientation dependence of /3 on {111} surfaces can 
be determined by first calculating the free energy F of a 
single step oriented at a fixed angle 0. To approximate 
this, we perform a low-temperature Ising expansion of the 
partition function, similar to the one used by Rottman 
and Wortisi^ They considered a step on a square lattice 
with one end fixed to the origin and the other end, a 
distance L away, fixed to the point {M = Lcos9, N = 
L sin 6) . Such a step is shown in Fig. QJl. The single-layer 
island (or compact adatom-filled region) is in the lower 
region, separated by the step edge — which is drawn as 
a bold solid line — from the upper part of the figure, 
representing the "plain" region. They found the step 
energies associated with the broken bonds of the step- 
edge to be 



E„ = e{M + N + 2n), n:^ 0,1,2,. 



(1) 



where e, sometimes^° called the "Ising parameter," is the 
bonding energy associated with the "severed half of the 
NN lattice-gas bond: Since the NN lattice-gas energy ei 
is attractive (negative) , and half of it is attributed to the 
atom on each end, it "costs" a positive energy e = — ^ei 
for each step-edge atom. Because longer steps require 
more step-edge atoms, the step energy is only a function 
of the step length: AI + N + 2n. Thus, Eq corresponds 
to the shortest possible step. To increase the length of 
this step, two more step-edge links - corresponding to 
one more step-edge atom - must be added, one going 
away from the fixed endpoint and one going toward it. 
Because this corresponds to two more broken bonds, in 
general, En+i — En ^ AE = 2e. With these energies, we 
can write down the partition function Zq, assuming 6 is 
fixed but L is large enough so that integer values of M 
and TV can be found: 



Ze = 5M,Jv(0)e-^°/*''^ +5M,Jv(l)£ 



-Ei/kT 



(2) 



where gM,N{n) corresponds to the number of ways a step 
of length M + N + 2n can be arranged between the two 
endpoints. 
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FIG. 1: There is a one-to-one correspondence between the 
shortest-distance steps connecting points on a square lattice 
[(a), left panel] and the shortest distance steps connecting 
points on a triangular lattice [(b), right panel]. This figure 
shows two corresponding steps. Analogous to the M + N 
broken bonds oriented at 0° and 90° on a square lattice, there 
are M' + N' broken bonds oriented at 0° and at 60° on a 
triangular lattice. However, there are another M' + N' broken 
bonds oriented at 120° on a triangular lattice. 



For low temperatures, only the first term in Eq. (j^J 
need be considered (Appendix I provides the lead- 
ing correction term, which gives a correction of order 
exp(— 2e/fcsT)). To lowest order, then, F is 



FfnEo- ksTln 



M + N 
M 



(3) 



where we have inserted the value of (7M,Ar(0) obtained 
from a simple combinatorial analysisi^iii After taking the 
thermodynamic limit {M, TV 3> 1) and using Stirling's 
approximation, F becomes 

FKEo-kBT[{M+N) ln(M + iV) - M In M - N In N] (4) 



B. Triangular Lattice Step Energy 

In order to extend this result to a step on a triangular 
lattice, we need to make a few minor adjustments. First, 
we introduce a linear operator C that transforms the co- 
ordinates of a point on a square lattice (M,7V) to those 
on a triangular lattice {M',N'); cf. Fig.^. This opera- 
tor finds the position of a point in a coordinate system 
whose positive y-axis is bent at 60° with respect to the 
positive x-axis: 



N' 



C 



1 -l/\/3\ (M 
2/V3 J \N 



(5) 



With the aid of £, we can see how Eq changes on a 
triangular lattice. To begin, we imagine a step in the first 
sextant (from 0° to 60° degrees in the plane) starting at 
(0,0) and ending at (M', TV'). Such a step is shown in 
Fig. nj). As before, the bold solid line represents the 
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step edge with the bottom region a layer higher than the 
top (or, alternatively phrased, it separates the upper, 
adatom-free region from the lower, adatom-filled region). 
The broken bonds required to form the step will have 
only three orientations: 0°, 60°, and 120°. If we consider 
the shortest step between the two points (corresponding 
to energy Eq), then there will be exactly M' + N' broken 
bonds oriented at 0° and 60° (these bonds are analogous 
to those oriented at 0° and 90° on a square lattice). There 
will be another M' + N' broken bonds oriented at 120° 
(drawn as bold, dashed lines in Fig. QJ. In total, there 
will be 2{M' + N') broken bonds. Since e is the energy 
of these severed bonds, = 2e(M' + N'). Thus, the 
energy is proportional to the step length, as was the case 
on a square lattice. Using £ to write M' and N' in terms 
of M and N gives 



2e M 



2sL cos 



sin 9 

71" 



(6) 



C. Triangular Lattice Step Degeneracy 

Next we consider the degeneracy factors g^{n) on a 
triangular lattice. For the ground state 5^(0) there is a 
one-to-one correspondence between the shortest distance 
steps connecting two points on a square lattice and the 
corresponding steps on a triangular lattice (see Fig. 
Therefore, we know that the degeneracy factor g^j, j^, (0) 
for steps of energy Eq on a triangular lattice must be 
identical to gM\N'{0) implicit in Eq. ©I More precisely, 
if we assume the point {M,N) is in the first quadrant, 
and {M' , N') is in the first sextant, then on a square lat- 
tice, shortest-distance step-links are oriented at either 0° 
or 90°, whereas on a triangular lattice they are oriented 
at either 0° or 60° (i.e. in the first sextant, the shortest 
path cannot have links oriented at 120°). In both cases, 
the individual step-links can only be oriented in one of 
two directions and, therefore, besides the transformation 
between coordinates, the total number of path arrange- 
ments is the same. 

Using Eq. Q and Stirling's approximation, we find 
the low-temperature free energy (Appendix I provides 
the lowest-order correction): 

F^E^~ keT ln[gM',JV'(0)] « E^ - fc^Tx 

X [{M' + N')hi{M' + N') - M'\tlM' - A^'lnA^'] (7) 

Alternatively, we can transform Eq. Q for the square 
lattice to the triangular lattice by just replacing N/M = 
tan6' with (2tan6')/(\/3 — tan6'). (This ratio is just 
N' /M' , so it might be termed tan 6*'.) [We must also 
make a simple (and ultimately inconsequential) change 
to ^0-] 

We now take the thermodynamic limit (M', A^' ^ 1) 
and write M' and A^' in terms of M = L cos 9 and A^ = 
L sin 9 via Eq. (O . Then dividing by L and definingiS 

sin 9 2 
r/± {9) = cos 6* ± , 7/0 {0) = ^ sin 6*, (8) 



all non-negative in the first sextant, we straightforwardly 
find the step-edge line tension (or free-energy per unit 
lengtbli) 13{9): 

a\\p{9) =. 2eT]+{9) - kBT[s+{9) - s^{9) - so{9)] , (9) 

where an is the nearest-neighbor spacing and 



s,{9) ^Tj,{9) In f],{9), z = +,0,-. 



(10) 



For the special case of the maximally kinked orientation, 
Eqs. © - (Cni) reduce to 



afc/3(30°) = 2£- fcsrin2. 



(11) 



where — (\/3/2)a|| for the {111} surface. This result 
for the maximally kinked case (including steps at 9 ~ 
45° on a square lattice) was derived earlier from a direct 
examination of entropyi^ 

For specificity, we recall some established results. For a 
hexagonal lattice with just nearest-neighbor attractions, 
the critical temperature Tc is long knownti^ 



ksTc = 2£/ln3 w 1.82£ 



(12) 



From the equilibrium shape of islands over a broad tem- 
perature range, Giesen et alm& deduced that the free en- 
ergy per lattice spacing in the maximally kinked direc- 
tions is 0.27 ± 0.03eV on Cu{lll} and slightly smaller, 
0.25 ± 0.03eV, on Ag{lll}. Combining these results 
with Eq. ini), we find e is 0.126eV on Cu{lll} and 
0.117eV on Ag{lll}. In both cases, then, room tem- 
perature is somewhere between Tc/9 and Tc/8. 



D. Main Result: Simple Expression for Low-T 
Stiffness 



As shown just above, the step-stiffness /3 = f3{9)+(5" {9) 
computed from Eq. Q depends to leading order only on 
the combinatoric entropy terms sq and s± of Eqs. 
and ((Tn|l . Hence, 



so thaU& 



Si — Si 'T S,j^ 



'2 2 



With this notation, the reduced stiffness is 



knT 



so + s- ~ s+, 



where 



n/3 ' 



V3 



So 



s± 



2 cos 26* 
71 sin6' ' 

-2cos26iT2\/3sin26' 
3 cos sm t) 



(13) 
(14) 

(15) 

(16) 
(17) 
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FIG. 2: (Color online) As the temperature drops close to 
Tc/Q (just below room temperature for Cu and Ag{lll} sur- 
faces), the numerical evaluation of the exact stiffness- ap- 
proaches the solid line representing the low-temperature ap- 
proximation given in Eq. IjlB^ . The small circles indicate eval- 
uations using the exact results of Eqs. 1231 and 1251 . At 
Tc/O, when 6 decreases, the exact solution begins to deviate 
from the approximation when its curvature changes sign near 
6 « 7r/100 = 1.8°. The scale here is linear, in contrast to the 
logarithmic scale of Fig. 2 of Ref. 0. The inset shows more 
fully how the exact stiffness approaches the low-temperature 
limit for the particular azimuthal angle S = 7r/24 = 7.5°. 



Adding these terms together gives our main result - a 
remarkably simple form for the reduced stiffness in the 
low-temperature (T <C e/ks) hmit: 



3™ — m 



/3a|| 

where m = t&nO. 



2x/3 2%/3(l +m2)3/2' 



(18) 



E. Synopsis of Exact Results and Application to 
Range of Breakdown of Low-T Limit Near ^=0 

To test how low the temperature should be for Eq. H18|) 
to be a good approximation, we compare it to a numer- 
ical evaluation of the exact implicit solution of the Ising 
model. The derivation of this solution, outlined by Zia^^ 
gives a 6*'* order equation for (3{9). In essence, after con- 
version to our notation, his key result for the step free 
energy /? is given bj*ii 



/3q|| 
knT 



^M0)MO,T/T,)+7j^{e)MO,T/T,), (19) 



where the ip^s are the solutions of the pair of simultaneous 
equations for the angular constraint. 



sinh?/;2 — si'nh(ipi — ip2) V- 
sinh?/;! + sinh('!/'i— ■!/'2) rjo 

and the thermal constraint, 



V3cot6' -1 



cosh ^fji + cosh ■02 + cosh(?/'i — ^p2) — f{z)'- 



l + 3z^ 

2(z--z2 



(20) 



(21) 



where z = exp{-2e/kBT) = ^-'^"/^ , via Eq. (O. The 
ratio r]-/r]o of Eq. (|20|l is a monotonically decreasing 
function which is oo at 9^0°, 1 at 61=30°, and at 9=60°. 

In these high- symmetry directions, Eqs. (|20|l and (|21|) 
yield analytic solutions for (3 and /?: 



keT 

/3(0)«|| 
ksT 



= 2cosh-r^i±^I±^) (22) 



2 V2(3 + 2/)(/-- V3 + 2/) 
3 



V3 



cosh 



V3T27-1 

-1 ^/-l' 



ksT 

/3(V6)a|t _ V3(/-3)(/-f 1) 



keT 



/ + 3 



(23) 
(24) 
(25) 



Details are provided in Appendix II. Akutsu and 
Akutsui^ also derived these equations, in different 
notationA2, and from the more formal perspective of the 
imaginary path-weight method. Symmetry dictates that 
the solution at 6* = 60° = 7r/3 is the same solution as 
that at 6* = 0°. Furthermore, at T = T^, f{z) = 3, so 
Eqs. (I22Jl-lE3) all go to 0, as expected. 

To find f3 in general directions, we solve Eqs. (|20|) 
and [or, equivalently, Eq. lIH^ ] numerically. As 

Fig. 121 shows, once T decreases to nearly Tc/9, Eq. ((TH|) 
more or less coincides with the exact numerical solution 
for the stiffness. At such low temperatures (compared 
to Tc), the approximation only fails below some very 
small, temperature-sensitive critical angle 9c- Although 
it might seem easy to determine this angle by eye, es- 




Tc/T 



FIG. 3: (Color online) Two estimates for the critical angle 6c, 
below which the approximation given in Eq. (1181 begins to 
fail, as a function of Tc/T. The black dots connected by the 
solid, blue line represent the first estimate, defined to be the 
angle 9i at which the curvature of the numerically determined 
inverse stiffness changes sign. The dashed, red line represents 
the second estimate, 02 = cot~^[(4/ — 1)/(50V3)]. At angles 
below 9c, the three theorems of Section UTTI break down, and 
higher order terms are required in the expansion of the step 
partition function. At temperatures between Tc/9 and Tc/8 
(roughly room temperature for Cu and Ag{lll} surfaces), Oc 
is on the order of a few degrees. 
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timating it quantitatively turns out to be a subtle and 
somewhat ambiguous task. We discuss two possible esti- 
mation techniques below. 

In the first approach, we estimate 9c to be the angle 9i 
at which the curvature of the exact solution changes sign. 
The points on the solid curve in Fig.|21show 6i at several 
temperatures ranging from Tc/9 to Tc/4:. At tempera- 
tures near and above Tc/4, 9i does not reliably estimate 
9c because there is a sizable curvature- independent differ- 
ence between the exact solution and the approximation 
given in Eq. ifT^ evident even a.t 9 = 30° (see Fig.EJ. On 
the other hand, as the temperature dips below Tc/5, this 
difference fades, and the use of 9i to estimate 9c becomes 
ever more precise. 

A second, more fundamental way to estimate 9c comes 
from an examination of the assumptions required to de- 
rive the simple expression for the low-T limit Eq. H18() 
directly from the exact solutions Eqs. H2()|l and (|21(l . In 
Appendix III we show that to do so 9 must be greater 
than some 92 satisfying 

cot02«(^^^). (26) 

To give definite meaning to this inequality, we estimate 9c 
directly from Fig.|21at a single temperature, say Tc/5. At 
that temperature, 9c is nearly 10°. If 6*2 is to accurately 
represent 9c, it should also be around 10° at Tc/5. We 
enforce this by interpreting the in Eq. (|26|l to mean 
'= 1/50.' The dashed (red online) curve in Fig.|21shows 
the resulting 02 as a function of temperature. Clearly 
9i and 92 are very different estimates for 9c- While 6*2 
is reliable at all temperatures (unlike 9i), it is less pre- 



^(30°)/;8(0°) >3(30°)/^(0°) 




Tc/T 

FIG. 4: (Color online) Ratio of the stiffness [solid red curve, 
left ordinate] and the free energy per length [dashed blue 
curve, right axis] for edges oriented in the maximally zig- 
zagged (6^ = 30°) and close-packed {9 = 0°) directions, based 
on taking the ratios of Eqs. 1251 and 1231 and of Eqs. 1241 and 
1221 1. respectively. The line-tension ratio increases slowly but 
monotonically to the r=0 limit 2/^/3 ~ 1.15. In contrast, 
the stiffness ratio plummets toward 0, the value predicted by 
Eg. 1)18^ . providing an indicator how low T must be for this 
simple low-T formula to be a good approximation at all an- 
gles. 



cise than 9i at lower temperatures. A combination of 9i 
and 92 is therefore the best estimate for 9c, being closer 
to 9i at lower temperatures, and closer to 92 at higher 
temperatures. 

In essence, 9c is no more than a few degrees between 
Tc/9 and Tc/8, regardless of which estimation technique 
is used. We therefore reach the practical conclusion that 
Eq. 1)18(1 is valid for almost all angles at temperatures 
near and below Tc/8, which fortunately happens to be 
around room temperature for Cu and Ag{lll}. 

Finally, we emphasize that f3 varies significantly with 
angle, especially at lower temperatures (where the equi- 
librium crystal shape (ECS) is hexagonal rather than cir- 
cular). If one wants to approximate (3 as isotropic rather 
than using Eq. (fT^ . one should not pick its value in the 
close-packed direction (viz. 9 = 0°); Fig. 0| provides stun- 
ning evidence of this conclusion. From Eq. H18() we also 
see that at low-temperatures the stiffness actually in- 
creases linearly with temperature. This contrasts with its 
behavior at high temperatures, where /3 must ultimately 
decrease as the ECS becomes more nearly circular and 
the steps fluctuate more easily. 



III. GENERAL RESULTS FOR STIFFNESS IN 
LATTICE-GAS MODELS IN 
LOW-TEMPERATURE APPROXIMATION 

In this section we present three theorems that are valid 
under two conditions: First, the energy term in the free 
energy must be a linear combination of cos 9 and sin 9. 
From Eq. © and [implicitly] Eq. (O we see that this 
property holds true in general for lattice-gas models, even 
when considering next-nearest neighbors and beyondiSi 
Second, the temperature must be low enough so that the 
entropy is adequately approximated by the contribution 
of the lowest order term, ks \ng{0). This entropic contri- 
bution is due exclusively to geometry or combinatorics of 
arranging the fixed number of kinks forced by azimuthal 
misorientation. Hence, it must vanish near close-packed 
directions (0° and 60° in the first sextant). For angles 
sufficiently close to these directions, in our case less than 
9c, the leading term becomes dominated by higher-order 
terms, and the three results no longer apply. 



A. No Contribution from Energy to Lowest-Order 
Stiffness (LOS) 

The first theorem is a remarkable consequence of the 
first condition, that the energy term in the free energy is 
a linear combination of cos 9 and sin0. Since the stiff- 
ness P{9) = f3{9) + (3" (9) and since cos" 6* = -cos 61 
and sin" 6' = — sin 6*, we see that the lattice-gas energy 
makes no contribution whatsoever to the low-T limit of 
reduced stiffness, as shown explicitly for square lattices 
long ago.^'^"'^ Thus, we retrieve the result that the leading 
term in a low-temperature expansion of the reduced stiff- 
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ness (3(9) /ksT depends only on g{0), which is determined 
solely by geometric (combinatoric) properties. Of course, 
higher-order terms (which are crucial near close-packed 
directions) do have weightings of the various configura- 
tions that depend on Boltzmann factors involving the 
characteristic lattice-gas energies. Furthermore, next- 
nearest-neighbor interactions can (at least partially) lift 
the (?(0)-fold degeneracy of the lowest energy paths.^ 



B. Step Line Tension Not Extractable from LOS 

An important corollary is that from the stiffness it is 
impossible to retrieve the energetic part of the step free 
energy, the major component of P{9) at lower tempera- 
tures when the islands are non-circular! Thus, contrary 
to a proposed method of data analysis, ^"'^ one cannot re- 
generate (3{9) from P{9) by fitting the stiffness to a sim- 
ple functional form and then integrating twice. In this 
framework, the linear coefficients of cos 9 and sin 9 can be 
viewed as the two integration constants associated with 
integrating a second-order differential equation^ 



C. LOS on fcc{lll} Has 6-fold Symmetry 

1. 1. General Argument 

Another important result is that the leading term in 
the stiffness at low temperature has the full symmetry of 
the 2D net of binding sites rather than the possibly lower 
symmetry associated with the full lattice. Specifically, 
for the present problem of the {111} face of an fee crys- 
tal, the stiffness (3{9) to lowest order has the full 6-fold 
symmetry of the top layer rather than the 3-fold symme- 
try due to symmetry breaking by the second layer. In 
contrast, the step energy of B-steps ({111} microfacets) 
differs from that of A-steps ({100} microfacets), leading 
to islands with the shape of equiangular hexagons with 
rounded corners, but with sides of alternating lengths 
(i.e., ABABAB). 

To see the origin of the 6-fold symmetry of the stiff- 
ness, suppose without loss of generality that steps in the 
X' direction have energy Ea per lattice spacing, so that 
those in the Y' direction have energy Eb- Furthermore, 
we must make the crucial assumption that any corner 
energy is negligible. Then all shortest paths to {M'.N') 
have the same energy M'Ea + N'Eb, with degeneracy 
still gM',N' (0). Thus, the free energy is M'Ea + N'Eb - 
kBTlngM' .N'{0), while that of its mirror point (through 
the line at 6* = 30°) is N'Ea + M'Eb - kBT\n gw' ,M'iO)- 
The crux of the proof is that gN',M'{0) = gM',N'{0)- 
Thus, while the free energies at the pair of mirror points 
differ, the energy parts are obliterated when the stiffness 
is computed (since M' and iV' are linear combinations 
of cos 9 and sin0), leaving just the contribution from the 
entropies, which are the same to lowest order. 



2. 2. Orientation- Dependent 3-Atom Interaction 

Within lattice-gas models with only pair interactions, 
there is no obvious way to distinguish A and B steps; 
the minimalist way to obtain different step energies for 
A and B steps within the lattice-gas model is to in- 
voke a non-pairwise 3-site "trio" interaction associated 
with three [occupied] sites forming an equilateral tri- 
angle with NN sides. In contrast to the ones consid- 
ered heretoforef^^'^** these novel trio interactions must be 
orientation- dependent: If the triangle points in one di- 
rection, say up, the interaction energy is positive, while 
if it points in the opposite direction, it has the opposite 
sign. (Of course, there could be a standard orientation- 
independent 3-site term in the Hamiltonian. As in the 
analogous situation for squares, we expect that such a 
term would simply shift the pair interactions, at least 
in the SOS approximation as discussed in Ref. 9.) The 
contributions from such a symmetry-breaking interaction 
would cancel in the interior of an island (in the 2D bulk), 
but would distinguish A and B edges. Specifically, each 
side of the equilateral triangle is associated with a link, so 
that 1/3 of its strength can be attributed to each. Each 
link has a triangle on both sides, one of each orientation. 
Hence, the difference between the energy per ay of A and 
B steps is 1/3 the difference between the trio interactions 
in the two opposite orientations. 

For the ground-state, minimum-number-of-links con- 
figurations, such a term will not lift the degeneracy since 
each configuration has the same 1) number of horizontal 
{X') hnks, 2) number of right-tilted diagonal finks {¥'), 
and 3) difference between the number of convex and con- 
cave "kinks" (i.e., bends). Since this statement is not 
true for higher-energy configurations, the 6-fold symme- 
try is not preserved at higher orders. Nonetheless, at low 
T it should be a decent approximation for the stiffness 
(much better than for the island shape). 

Thus, our result that the breaking of 6-fold symmetry 
on an fee {111} is much smaller for the stiffness than for 
the free energy, is more general than the nearest-neighbor 
lattice gas model which underlies Eqs. © - (|10|l and the 
resulting Eq. H18I) derived below. We reemphasize that 
the necessary assumptions are 1) that the orientational 
dependence of the step energy be just a linear combina- 
tion of sm9 and cos 9 and 2) that no interaction break 
the degeneracy of the shortest path corresponding to ori- 
entation 9. As above, for angles near close-packed direc- 
tions, the higher-order terms become important at lower 
temperatures than for general directions. This feature is 
illustrated in Fig. Inland its associated formalism is given 
in Appendices I and III. 



IV. COMPARISON TO EXPERIMENT 

In Fig. |31 we compare Eq. H18|) to measurements on 
Cu{lll} and Ag{lll}. The experimental data were de- 
rived from the equilibrium shape of 2D islands using the 
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FIG. 5: (Color online) A comparison between Eq. 1181 and 
experiments on Cu and Ag{lll}. Eq. 11811 appears as a solid 
black line, while the average of the experimental data is a 
thick dashed blue line. The agreement is reasonable at all 
angles. In either case the thin dashed red line is a [smoothed] 
average of the data for the given angle while the thin solid 
purple line corresponds to the angle mirror-reflected through 
a radial at 30°, i.e. at 60° — 0. 



method described in Ref. IQ. The solid black line corre- 
sponds to Eq. lfTH|l . while the thick dashed blue line corre- 
sponds to the average of the experimental measurements. 
Eq. (|18|l captures the overall trend and is satisfactory at 
most angles and temperatures. As expected, Eq. |[TH|) 
somewhat overestimates (3 near 6* = 0° (since the T — Q 
singularity remains) . Furthermore, near 6 = 30° Eq. H18() 
somewhat underestimates the experimental $, but only 
by a factor of 1/6 for Cu{lll} and 1/4 for Ag{lll}. 
This is in striking contrast to the analogous NN theory 
for Cu{001} near 45°, which underestimates /3 by a fac- 
tor of 4. Finally, notice there is no clear temperature 



dependence in the measured data. This is further evi- 
dence that fS/ksT is a constant at low-temperatures, as 
Eq. l(TH|) suggests. 

The agreement between theory and experiment is 
a pleasant surprise considering analogous comparisons 
made for Cu{001}- found /3 to be four-times larger than 
the theoretical value at large angles (near 9 = 45°). It 
was later shown^iM, that this discrepancy could be par- 
tially accounted for by considering next-nearest-neighbor 
(NNN) interactions (or right-triangle trio interactions, 
which turn out to affect (3 at low temperatures in the 
same way) . Clearly, the success of Eq. lfTH|l suggests that 
these interactions are less relevant for {111} surfaces. 
This is reasonable because the ratio of NNN distance to 
NN distance is smaller by a factor of a/2/3 on a triangu- 
lar lattice compared to on a square lattice. Furthermore, 
in the close-packed direction {9 ~ 0°), for every bro- 
ken NN-bond there are only one and a half broken NNN 
bonds on a triangular lattice, compared to two broken 
NNN bonds on a square lattice. These simple arguments 
help explain why NNN-interactions may increase (3 by 
only 20 to 30% on Cu/Ag{lll}, as opposed to 400% on 
Cu{001} surfaces. 



V. CONCLUDING DISCUSSION 

By generalizing the low-temperature expansion of the 
nearest-neighbor square lattice-gas (Ising) model to a tri- 
angular lattice, we have found a remarkably simple for- 
mula for the orientation dependence of the {111} surface 
step stiffness. This formula, unlike its square lattice ana- 
log, fits experimental data well at general angles, suggest- 
ing that NNN-interactions are relatively unimportant on 
{111} surfaces. 

To corroborate this picture and explain the success of 
Eq. (|18() . we are currently using the VASP packaga^to 
perform first-principle calculations. In particular, we 
are examining the ratio of the NNN to NN interaction 
strength. Preliminary results^^ suggest that this ratio 
is roughly an order of magnitude smaller on Cu{lll} 
than on Cu{001}, and essentially indistinguishable from 
zero. This tentative finding is consistent with expecta- 
tions from the semiempirical embedded atom method, 
which predicts that indirect interactions are insignifi- 
cant/negligible between atoms sharing no common sub- 
strate atomsi^ We are also calculating the difference 
in trio interactions between oppositely oriented triangle 
configurations. 

We expect that our formula, as well as the general 
6-fold symmetry of the stiffness (except in close-packed 
directions), should be broadly applicable to systems in 
which multisite or corner energies are small and for which 
the bond energies are considerably higher than the mea- 
surement temperature. Studies which ignore the 3-fold 
symmetry breaking on metallic fee {111} substrates, such 
as a recent investigation of nanoisland fluctuations on 
Pt{lll})2i should be good representations. Many recent 
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investigationg^"^^ focus on the larger asymmetry of the 
kinetic coefficient^^ taking the stiffness to be isotropic. 
In such cases, this stiffness should not be characterized 
by its value in close-packed directions. 



APPENDIX I: LEADING TERM IN 
LOW-TEMPERATURE EXPANSION 

A. Review of Results for Square Lattice 

In this appendix we discuss the lowest-order correction 
to the ground state entropy of the step running from the 
origin to an arbitrary particular point. First we review 
results for a square lattice. We can rewrite Eq. (|2Jl as 



^ ^ 9M,n{1) ^-AE/kT 

9mm (0) 



(27) 



Then, assuming the exponential is small, we have 
F^Eo- ksT (ln[5M,^(0)] + iM^e-'^/>^sT ^ ^^^8) 

L 5M,Ar(0) 

A combinatorial analysis^iii*^^ shows that 

9m, nI^) 
Then Eq. Q generalizes to 



F w Eo- kBT[{M + N)\n{M+N)- M\nM - N\nN 



-2e/kBT'_ 



MN 



(30) 



B. Results for Triangular Lattice 

For the triangular lattice we find important differences 
from the square lattice for the higher-order terms. Specif- 
ically, we consider how ^(l) changes. In contrast to 5(0), 
we cannot simply replace M and N with Af and N' . 
There is no one-to-one correspondence between paths 
of energy Ei on a square lattice and those of energy 
E^ on a triangular lattice. This failed correspondence 
for higher terms follows from the observation that E^- 
configurations are only one link longer than iSp -steps, 
whereas £'i-configurations are two links longer than Eq- 
steps: - E^^ = AE^ = e, or 



2M 



n = 0,1,2, 



(31) 



Hence, we require a separate combinatorial analysis. 

We imagine a step of energy Ei in the first sextant. 
Such a step (see Fig. EJ will have either: (1) (Af -|- 1) 
links oriented at 0° (denoted "X-hnks"), (A^' - 1) hnks 
oriented at 60° (denoted "Y- links" ) , and one link oriented 




(M',N') 



FIG. 6: Two equivalent steps having energy Ei . The dashed 
step contains (Af' -I- 1) X-links, (A*'' — 1) Y-links, and one T- 
link, while the solid step contains {M' — 1) X-links, (A'^' + 1) 
Y-links, and one T-link. 



at 120° (denoted "T-links"), or (2) (M'-l) hnks oriented 
at 0°, (A^' + l) links oriented at 60°, and one link oriented 
at —60°. In the first case, the problem can be reworded 
as follows: how many ways to arrange an {M' + N' + 1)- 
lettered word with (Af + 1) X's, {N' - 1) Y's, and one T. 
In the second case, the problem is the same, only with 
M and N switched. Thus, the solution of this traditional 
combinatorial problem gives the total number of next-to- 
shortest paths (7^(1): 



M'+N' + l 
M' + l 



M'+N' + l 
N' + l 



M'. (32) 



With 5'^ (0/1) and E^ in hand, we can write the low- 
temperature partition function expansion for a triangular 
lattice. Using Eq. © and expanding the logarithm as in 
Eq. ((2H1), we have 



F 



E^ - ksT \ ln[.g(0)] + ^e-^^^^.T 



(33) 



Taking the thermodynamic limit (Af, TV' 1) and using 
Stirling's approximation gives 

F « E(^ ~kBT[{M' +N')\tl{M' + N')-M' hiM' -N' \nN' 
./J, rr.M'^ + N'^ + M'N''^ + N'M'^ , , 



-I- e~ 



M'N' 



The pair of cross-factors in the last coefficient are absent 
in Eq. (|Sn|l for the square lattice. 

The correction term becomes non-negligible when the 
final term in Eq. H34|l becomes of order unity. At low T 
this occurs only near close-packed directions, so for small 
values of 9. In this regime, to lowest order in 6*, N' = 
(2L/V3) sin e ^LQI-H and M' = LcosO- N'/2 -> L. 
Then the critical value of 9 is 



Specifically, based on Eq. P5|l and using e R 

Cu{lll}, we find that 9i^^ is 0.353°, 3.18 
for T/Tc of 1/9, 1/5, and 1/4, respectively. As clear 



(35) 

0.12 eV for 
, and 5.51° 
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from Fig.O this criterion turns out to underestimate the 
values for 6c obtained in Section II. D, mainly because 
Eq. H^^5(l was derived from an expression for f3{6) instead 
of f3{6) (which should depend more sensitively on i 
the plotted thermal range. 



over 



APPENDIX II: EXACT FORMULAS FOR LINE 
TENSION AND STIFFNESS IN MIRROR 
DIRECTIONS 

A. General results for all orientations 

In this appendix, we derive Eqs. (|22|l - (|25|l for the 
mirror-line directions 9 — 0° and 9 — 30° from Zia's 
implicit exact solutioni^ 

To begin, because /3 =^ (3 + (3" (where the prime rep- 
resents differentiation with respect to 9), it follows from 
Eq. ^ that 



(36) 



We can simplify Eq. (|36|l by finding relationships between 
the various derivatives of the V''s. Differentiating Eq. H21(l 
with respect to 0, regrouping, and using Eq. (|20|1 . we get 



+ -02 '7- 
Differentiating again yields 



(37) 



(38) 



Using Eq. H38|) , we rewrite the last part of Eq. H3t)|l (con- 
taining -01 and V'2) in terms of just ip'^ and ■02 ■ Then, 
using Eq. H37() we eliminate 02 in favor of -01 • We are 
then left with an equation relating /3 to only 0^: 



0i 



201 



\/3 cos 9 — sin 9 



(39) 



For general angle, we must evaluate tpi numerically. 
However, for the two high-symmetry directions we can 
obtain analytic results that allow us (with the aid of 
Eq. (|19|l for f3) to write explicit expressions for P, as 
presented in the next two subsections. 



B. Results for 9^0° 



At 9 = 0°, Eq. (Uni) reduces to 
Pa\\ 



02(0), 



(40) 



assuming that '01 (0) is finite. Furthermore, near 9 = 0°, 
Eq. H20|) can be inverted and the sinh's combined to get 



sinh(-0i - 1-02) cosh(i02) 2 
sinh('02 - 501 ) cosh(ii/ii) ^3 



(41) 



For Eq. (gll) to hold at 6* = 0°, 

02(0) = 2V^i(0). 
Eq. 12 1|) therefore becomes: 

2 cosh 01 (0) -f cosh(20i (0)) 



(42) 



(43) 



Solving this for cosh -01(0) and taking the positive root, 
we find: 

coshV'i(O) = cosh(iV;2(0)) = i(-l + V3T27), (44) 

consistent with the assumption of finite 0i(O). Solving 
for 02(0) and combining with Ec^. (|in|l yields Eq. (|^ . 
Correspondingly for (3, aX 9 = 0° Eq. H39|l becomes 



while Eq. (|S7|l becomes 

02(0) = 0, 



(45) 



(46) 



provided 01(0) is finite. We obtain 01(0) by differentiat- 
ing Eq. H41|) with respect to 9 and then setting 9 = 0° so 
that Eqs. 1021 and apply. This give 

-01 (0) = ^ tanh 0^1 (0) [1 + 2 cosh 0i (0)] . (47) 

By combining this with Eq. (|44|l for cosh?/;i(0), we see 
that 01(0) is indeed finite, as we earlier assumed. Thus, 
Eq. becomes Eq. if^ . as desired. 



At 9 



C. Results for 6 = 30° 

30° = 7r/6, Eq. (0 becomes 

/3(V6)a|| 



knT 



^[0i(V6)+0^('^/6)]. 



(48) 




FIG. 7: (Color online) Numerical evaluation of %p\ (dashed 
red curve) and 1/12 (solid blue curve) as functions of angles 
at temperature Tc/S equivalent to room temperature for the 
experimental systems Cu and Ag {HI}. Note that the linear 
behavior near one limit and the divergent slope near value 
zero at the other. At higher temperatures the curves are 
qualitatively similar but progressively smaller in magnitude. 
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Furthermore, near 9 = 7r/6, rjo/r]^ « 1 + 2\/3 A9, where 
A9 = 9 — 7t/6. Inverting Eq. (|20|l we therefore have 



sinhipi + sinh('i/;i — 1P2) 
sinh7/)2 — sinh(V'i — ^/'2) 



I + 2V3 A9, 



(49) 



By inspection, at 9 = tt/G {A9 = 0), one solution to this 
equation is just 

MV6) = M''/6)- (50) 

Plugging this result into Eq. (|21(l and solving for ■ip2{'^/6) 
gives. 



cosh V'2 (Ve) 



/-I 



(51) 



Combining this with Eq. H48II (where we now know 
iji{n/6) = ■02(7r/6)) results in Eq. if^ . 
As for P, at 9 = tt/Q, Eq. ((23 becomes 



while Eq. (|37|l becomes 



2^-1 (Ve), 



-V^sWe). 



(52) 



(53) 



Like before, we can find ^P'i{t:/6) by differentiating 
Eq. H49|) with respect to 9. Taking the result and set- 
ting 9 — 7r/6, so that Eqs. (|50|) and H53() apply, gives 



\/3 sinh ipi 



(54) 



cosh -01 + 2 

Finally, we combine this result with Eq. H51(l for 
cosh 02(^6) = cosh0i(7r/6) andEq. to get Eq. 
as desired. 



Solving this pair of equations for e'^'^ and e'''^ gives 



^ 2/(z)?7o 
Vo+V-' 



2f{z)r,- 
Vo + V- 



(58) 



If we then replace f{z) by its low-temperature limit, 
l/(2z), Eq. lUni) becomes 



/3Q|| 



= In 



r/_ In 



V- 



(59) 



By noting rjQ + rj^ = ^+1 f-nd using the definition for z, 
Eq. (|59|l can be easily simplified to Eq. @, from which 
Eq. H18|l for /? was derived. 

By deriving the approximation given in Eq. @ (and 
thus Eq. (|18|l 'l in this way, we can determine when the 
approximation becomes invalid. Specifically, we require 
cosh 01^2 ^1- As we showed, the more restrictive 
of these inequalities is the one involving cosh^i, since 
cosh ipi is necessarily smaller than cosh 02 in the first sex- 
tant by a factor of 770 /rj- (which is less than 1). Thus, the 
main assumption is cosh -01 ^ 1, which, from Eq. I|58|l . 
is just 



2/(z)r?o 

Vo + V- 



> 1. 



(60) 



The solution to this equation, which we call ^2, is given 
by the following inequality: 



APPENDIX III: REDERIVATION OF EQ. 
FROM EXACT SOLUTION 

In this appendix, we re-derive Eq. (|18|l directly from 
the exact solution for f]{9) given in Eqs. H20() and (|21|l . 
To do so, we just assume cosh 02 S> il-fvo (remember 
that r]-/T]o decreases from 00 at = 0° to 1 at 6* = 30°, 
so that, between these angles, this condition also implies 
that cosh -02 ^1)- In this case, Eq. (|20|l can be solved 
to give 



cosh ■01 sa — cosh ■02- 

V- 



(55) 



Thus, if cosh -02 ^ V-fVo > 1: then cosh^Ai 3> 1. We 
show here that these assumptions for cosh -01 2, together 
with the low-temperature replacement of f{z) by l/(2z) 
in Eq. (|2U, are enough to derive Eq. (fT5|l . 

When cosh0i,2 ^ 1, then cosh?/)i,2 ~ sinh?/;i,2 ~ 
e'/'i>2/2. With these approximations, Eqs. ^ and (EH) 
become remarkably simple; 



(56) 
(57) 



cot 6*2 < 



4/-1 
V3 ■ 



(61) 



Because cot 9 decreases from 00 at 61 = to l/-\/3 at 
9 = 7r/6, we know that angles in the first sextant that are 
greater than 6*2 will also satisfy the inequality in Eq. H61|l . 
Thus, Eqs. © and H18|) are valid in the first sextant at 
all angles above 6*2 ■ 
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